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We study the half-filled Hubbard model on the geometrically frustrated face centered cubic (FCC) 
lattice, using an auxiliary field based real space technique. The low temperature state is a para- 
magnetic metal at weak interaction, an antiferromagnetic insulator (AFI) with fiux like order at 
intermediate interaction, and an AFI with 'C type' order at very strong interaction. Remarkably, 
there is a narrow window between the paramagnetic metal and the AFI where the system exhibits 
spin glass behaviour arising from the presence of disordered but 'frozen' local moments. The spin 
glass state is metallic at weaker interaction but shows crossover to pseudogap behaviour and an 
insulating resistivity with growing interaction. We compare our results to available experiments on 
FCC and pyrochlore based materials and suggest that several of these features are typical of three 
dimensional correlated systems with geometric frustration. 



The presence of geometric frustration in an interact- 
ing electron system disfavours long range magnetic order 
and promo tes a complex electronic state with short range 
correlation^^. Such effects have been intensely explored 
in t he c ontext of the quasi two dimensional (2D) organic 
salt^^ where, in some compounds, the triangular lattice 
structure gives rise to a spin liquid rather than conven- 
tional Neel order. 

Unlike in 2D, there is no organised body of work 
probing the interplay of geometric frustration and Mott 
physics in three dimensions (3D). There are intriguing ex- 
perimental results on disparate system^^HsU^ whose com- 
mon features do not seem to have been noticed. The 
3D frustrated Mott systems are realised on face centered 
cubic (FCC) and pyrochlore lattices. They both involve 
corner sharing tetrahedra, disfavouring simple Neel or- 
der in the insulating phase. The FCC examples include 
the cluster compoundJ^"^ like GaTa4Se8, GaNb4Se8, etc, 
some alkali fuUeride's of the form^ A3C60, and the 'B 
site ordered' double perovskited^^t^^, e.g., Sr2lnRe06. 
The pyrochlore examples include the molybdateJ^^^^^ 
R2M02O7 and iridatePH2llLn2lr2 07. Most of these ma- 
terials, at ambient ' pressure ^^Si, are insulators close to a 
Mott transition. 

While there is great variation among these materials, 
the following features seem to be shared: (a) . In the Mott 
phase they usually exhibit no long range order down to 
the low est tempe rature, sometimes with a hint of 'spin 
freezing'^ISEIinHlll, (b). On pressure driven metallisa- 
tion, the resistivity is very large but finite'' ^'^ at low 
temperature over a wide pressure window, and exhibits 
a negative temperature derivative, befo re ev entual 'nor- 
mal' behaviour, (c) Optical conductivitjl^'^il^ where avail- 
able across the pressure driven transition, indicates large 
transfer of spectral weight, (d) The Hall conductance^'^ 
has a spontaneous anomalous contribution indicating sig- 
nificant non coplanar character in the magnetic back- 
ground. In addition, (e). some of th ese sys tems exhibit 
superconductivity at low temperatur e l ^"^^ !. 

The complexity of the three dimensional structures has 
prevented a detailed exploration of the possible magnetic 



phases and the Mott transition in these systems. While 
each material involves its specific electronic model, in 
this paper we focus on the single band Hubbard model 
at half-filling on the FCC lattice. We use an approach 
that captures the unrestricted Hartee-Fock (UHF) state 
at zero temperature, but retains the crucial thermal fluc- 
tuations of the interaction induced 'local moments', and 
their impact on electronic properties. Using a combina- 
tion of Monte Carlo (MC) and variational schemes we 
establish the following: (z) In the ground state, increas- 
ing interaction leads, successively, to transition from a 




FIG. 1: The U ~ T phase diagram of the Hubbard model 
on the FCC lattice at half-filling. The ground state can be a 
paramagnetic metal (PM), with no local moments, a spin glass 
metal (SGM) or spin glass insulator (SGI) with disordered 
local moments, and an antiferromagnetic insulator (AFI) . The 
AFI has 'flux' like order upto U/t « 30, beyond which it has 
'C type' order. The classification into 'metal' and 'insulator' is 
based on the temperature dependence of the resistivity, p{T). 
dp/dT > indicates a metal, dp/dT < an insulator. At 
finite temperature the system also has a paramagnetic (Mott) 
insulating (PI) phase. The magnetic transition temperature, 
Tc, and the spin glass freezing temperature, Tsa (see text), 
are indicated. We show the extrapolation of the Tc ~ 0.6t^ /U 
asymptote, that describes the C//t 3> 1 transition, to highlight 
the large deviation from the short range Heisenberg result. 
The PG region involves a pseudogap in the density of states. 
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paramagnetic metal (PM) to a spin glass metal (SGM), 
a spin glass 'insulator' (SGI), and then an antiferromag- 
netic insulator (AFI). The AFI has flux like order at 
weaker coupling, and 'C type' order (see Supplement) 
in the very strong coupling Heisenberg limit, (m) The 
Tc scales in the ordered phase, as well as the notional 
glass transition temperature Tqq, are a tiny fraction, ^ 
a few percent, of the hopping scale due to the frustration. 
{iii) While the AFI phase has a clear gap and divergent 
resistivity as temperature T — > 0, the spin glass insula- 
tor has a pseudogap (PG) in the single particle density 
of states (DOS), non Drude optical response, and large 
hut finite resistivity at T — Q with dp/dT < 0. (iv) The 
transport trends match remarkably with experiments on 
FCC Mott systems, and allow us to make predictions 
about their magnetic and spectral properties. 

There is surprisingly little theoretical work on the mag- 
netic phases or the Mott transition in the FCC lattice. 
There is a very early calculation exploring a restricted set 
of mean field statea^Zlbut, in contrast to two dimensions, 
there does not seem to be any cellular dynamical mean 
field theory (C-DMFT) result handling the combination 
of correlation and frustration. We use an approach, sug- 
gested long back by Hubbard himself^*, for the model: 

{ij)a- I I 

The tij = —t for nearest neighbour hopping on the FCC 
lattice. We will set t = 1 as the reference energy scale, /i 
controls the electron density, which we maintain at n = 1. 
J7 > is the Hubbard repulsion. 

We use a Hubbard-Stratonovich (HS) transformation^^ 
that introduces a vector field mi(T) and a scalar field 
4>i{T) at each site to decouple the interaction. This 
decompositiorPsHni retains the rotation invariance of the 
Hubbard model and reproduces UHF theory at saddle 
point. We treat the vHi and (pi as classical fields, i.e, 
neglect their time dependence, but completely retain the 
thermal fluctuations in m.i. is treated at the saddle 
point level, i.e, <j), (0^) = {U/2){{n,)) = U/2 at half- 
fllling, since charge fluctuations would be penalised at 
temperatures T U. Retaining the spatial fluctuations 
of nii allows us to estimate Tc scales, and access the cru- 
cial thermal effects on transport. We will discuss the 
limitations of the method later in the paper. 

With this approximation the half-filled Hubbard prob- 
lem is mapped on to electrons coupled to field, nii. 

Heff = Yl *y 4<TCj<T-AA^-^ XI ^i-^i + j X (2) 

where rhu = fi — U/2. We can write this as i?e// 
H,i{ni,} + Hcu where H^i = {U/^)Y.^'^^■ The 
{mi} configurations follow the distribution P{m.i} (x 

Within the static HS approximation H^ff and P{mj} 
define a coupled fermion-local moment problem. This is 



similar to the 'double exchange' problem, with the cru- 
cial difference that the moments are self generated (and 
drive the Mott transition) rather than fixed in size. Due 
to the fermion trace, P{m,i} is not analytically calcu- 
lable beyond weak coupling. To generate the equilib- 
rium {m.i} we use Monte-Carlo sampling. Computing 
the energy cost of an attempted update requires diago- 
nalising Hei- To access large sizes within limited time, 
we use a cluster algorithm for estimating the update cost. 
We calculate the energy cost of an update by diagonaliz- 
ing a cluster (of size Nc, say) around the reference site. 
We have extensively benchmarked this 'traveling clus- 
ter' methocP^. The static HS approach, retaining spatial 
fluctuations, has found successful application in corre- 
lated systems befor^2ll, xhe MC was done for lattices 
of size upto N = 12 x 12 x 12, with clusters of size 
Nc = 4x4x4. We calculate the thermally averaged 
structure factor S{q) = ' mj)e*'i'^''*~''j^ at 

each temperature. The onset of rapid growth in iS'(q) 
at some q = Q, say, with lowering T, indicates a mag- 
netic transition. Electronic properties (see Supplement) 
are calculated by diagonalising Hd on the full lattice for 
equilibrium {m.^} conflgurations. Since the MC ground 
state can be affected by annealing protocol, wherever pos- 
sible we have tested it against variational choices of {m^}. 

Fig{l] shows the U — T phase diagram of our model. 
First focus on the magnetism at T = 0. (i) The MC based 
minimization, j^(i?e//{nii}) = 0, leads to a state with 
rui = |mi| = for U < Uci ^ -it. (ii) For Uci < U < Uc2, 
where Uc2 ~ 6.7i, the ground state involves flnite m^, 
with a flnite width distribution P{m, U) = {j^ J2i '^("^ ~ 
I m.i I)), but with no long range spatial correlation. The 
system behaves like a spin glass with short range 'flux 
like' correlations. (Hi) Beyond Uc2 the ground state has 
long range flux like order till U ^ 3Qt, beyond which the 
virtual hopping generated exchange is effectively nearest 
neighbour and we obtain 'C type' order. The C type 
order is indeed expectecP^ for the AF Heisenberg model 
on the FCC lattice. 

The thermal physics deep in the AFI phase is con- 
trolled by angular fluctuations of the local moments, m^, 
about the ordered state. For U Uc2 this leads to 
the usual Tc oc t'^/U, but with a coefficient of w 0.6, 
much smaller than ~ 1.4 in the simple cubic case. At 
weaker interaction, for Uc2 < U ^ 30t, longer range 
and multi-spin couplings between the m^ become rel- 
evant and the Tc deviates significantly from the t'^/U 
asymptote. For U ^ 9t the size of the local moment it- 
self diminishes rapidly, due to increase in itinerancy, and 
the Tc falls sharply. Below Uc2 ~ 6.7< where we have 
a glassy phase we make a crude estimate of the 'freez- 
ing temperature' from the MC based local relaxation 
tim^P, tUT,U) = (l/7V)Ei/o'"°^di(m,(0).m,(i)). If 
the system undergoes an ordering transition, on lower- 
ing T, there is a rapid growth in To„ accompanied by a 
growth in the structure factor S{q) at the q's associated 
with long range order (LRO). For a glass transition, one 
observes similar growth in Tq.,, , without any signatures in 
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FIG. 2: Colour online: (a) Temperature and U dependence 
of the resistivity. The U/t values are marked in the plot. In 
the PM window, U < Ud, the zero temperature resistivity 
p(0) vanishes and dp/dT > at low T. In the SGM phase, 
Uai <U < Uc, p(0) is finite, with dp/dT > 0. In the SGI 
phase, i.e, Uc < U < Uc2, p(0) is finite, rapidly grows with U, 
and shows dp/dT < 0. For U > Uc2 where the ground state 
has long range order and a gapped spectrum, p(0) is infinite, 
and dp/dT < 0. For the weakly insulating ground states, {U — 
Uc)/Uc 1, increasing temperature leads to a crossover to 
dp/dT > beyond a temperature we call Timt{U). (b) The 
variation of the p(0) with U /t. (c) The variation of the average 
moment rriav at T = with U. Our resistivity is measured 
on the scale of po ~ hao /ttc^ , where ao is the lattice spacing. 
For ao ~ 3A it will be ~ GOpflcm. 

S{q). For U > Uc2 we observe LRO as well as a rapid 
increase in r^y at a single temperature Tc(U). For the 
window Uci < U < Uc2, however, Tav rises, at a temper- 
ature we call Tsg{U), without associated LRO. We have 
also 'heated' the system up from T = and discovered 
that any assumed ordered state is quickly destabilized 
while the moments themselves survive. Tsg varies in the 
manner shown in Figjl] vanishing for U < Ud where 
there are no local moments. 

Figj2] shows the resistivity p{T,U). For U < Ud the 
resistivity p{T = 0) = 0, and dp/dT > 0. For U > 
Uc2 the system is gapped at T = 0, p(T = 0) — > oo 
and dp/dT < 0. These are the obvious metallic and 
insulating behaviour that one expects across a correlation 
driven transition. For Ud < U < Uc2, however, the 
T = resistivity is finite, with dp/dT > for Ud < 
U < Uc, where Uc ~ 5.8t, and dp/dT < for [/^ < [/ < 
Uc2- This behaviour would usually not be expected in 
a translation invariant system, and arises because of the 
scattering of electrons from the 'frozen' local moments. 
The growth of mav{U), the mean magnitude of rn^, leads 
to the enhanced scattering with increasing U and finally 
the divergence of p{Q) due to the opening of a gap. The 
variation of p(0) and ruav with U/t are shown in panels 
(b) and (c) respectively. We characterize the system as 
metallic, at a given T and [/, when dp/dT > 0, and 
insulating when dp/dT < 0. With this convention, an 



'insulator' may have finite spectral weight at w = in 
the optical conductivity o'{uj). 

While the U > Uc2 system would have a gapped DOS, 
and the U < Ud case is likely to have a featureless DOS, 
the glassy window in between may have unusual spectral 
features. Fig(3|(a)-(b) shows the DOS for varying U/t a,t 
T = 0.03t and T = 0.20t, respectively. For lack of space 
the colour codes for U/t are marked in panel (b) only. 
In panel (a), for U < Ud the DOS is featureless, but for 
5.4 < U/t < 7.5 it displays a PG, and for higher U/t 
there is a clear gap. The large U/t phase is magnetically 
ordered at this temperature. At the higher temperature 
in panel (b), where there is no trace of magnetic order, 
the PG feature extends over a much larger U/t window. 
The weaker U 'metals' in (b) have a deeper PG compared 
to panel (a), while the weak gap insulators now have a 
PG feature rather than a hard gap. The evolution in 
panel (b) essentially illustrates the paramagnetic Mott 
transition on the FCC lattice. 

Panels |3](c)-(e) show the T dependence of the DOS 
for three typical U/t in the glassy window, where the 
ground state has frozen local moments. They all share 
the feature of thermally induced loss of low frequency 
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FIG. 3; Colour online: Density of states N{ijj). Panel (a) 
shows the DOS at T = 0.03i for increasing U /t showing the 
crossover from a correlated metal to an AFI (with weak sur- 
viving order) through a wide pseudogap window. The colour 
code for U/t is marked in panel (b). Panel (b) shows the DOS 
at r = 0.20t where the crossover is between the PM and a 
PI through a much wider pseudogap window. Panels (c)-(e) 
show the temperature dependence at three fixed U in the spin 
glass window. Notice the emergence of a thermally induced 
PG at weak interaction, U — 5.2t, and the presence of the 
PG at T = itself for U = 6.2t and 6.7t. Modest changes 
of temperature, T ~ O.lt, leads to large asymmetric shift of 
spectral weight from a; ~ to tj ~ [/. 
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FIG. 4: Colour online: Optical conductivity (t(w). Panel (a) 
shows (j{uu) at T = O.OSt for increasing U/t (marked in the 
plot). The response evolves from a Drude character to the 
gapped spectrum of the AFI through a non Drude regime. 
Panel (b) shows similar evolution at T = 0.20t, from the PM 
to a PI. There is no Drude peak visible down to U = 5t, the 
peak locations have moved to higher oj, and the overall scale 
of (j{uj) is halved. Panels (c)-(e) show the temperature de- 
pendence of (j{iu) for some U/t in the spin glass window. The 
weak moment system in (c) shows essentially a broadening 
Drude response with increasing T. The larger moment system 
with U/t — 6.2 in (d) shows a non Drude response with low 
frequency weight suppressed by temperature for T > 0.03t. 
Panel (e) shows a spin glass with large T = resistivity. The 
very low frequency weight (on the scale of the pseudogap) 
increases with T, while the weight at a; ~ [/ reduces with 
increasing T. The scale ao = 1/po- 



weight which shows up at a; ^ [/. There is markedly 
less change with T on the negative frequency side, par- 
ticularly in panels (d) and (e), compared to positive fre- 
quencies. This is due to the large asymmetry in the tight 
binding DOS of the FCC lattice. There is a subtle low 
energy difference between panels (c)-(d) and panel (e). In 
(c)-(d) the loss in low frequency weight, within w ^ ±i is 
monotonic with T. In panel (e), however, which neigh- 
bours the AFI, the low frequency weight first increases 
with T, upto T ^ O.lt, and then again diminishes at 
higher temperature. 

Fig|4] shows the optics for the same parameter choice 
as the DOS plots. Panels (a)-(b) show the evolution of 
(t(cli) across the metal-insulator transition, between the 
PM and AFI at T = 0.03t, and between the PM and PI 
at T = 0.20i. There is a clear window of non Drude re- 
sponse at low T, roughly corresponding to the PG regime 
m Fig|3l(a). In|4](b) the non Drude window in U/t has 
increased as in Fig|3](b) with a general suppression in 
the magnitude of cr(u}). The panels (c)-(e) show the sup- 



pression of low frequency optical weight, with some of 
it appearing at w > C/. Unlike the single particle DOS, 
the total optical weight is not conserved and varies with 
the kinetic energy. At U/t = 6.7, the very low frequency 
optical response is non monotonic in T, showing a quick 
increase and then a gradual suppression. This directly 
relates to the behaviour of p{T) in Figj2](a). 

We have highlighted a host of magnetic, transport and 
spectral features associated with Mott phenomena on the 
FCC lattice. However, like all many body methods, our 
approach too is approximate and let us touch upon the 
possible shortcomings before we atte mpt to relate our 
results to experiments. Earlier paperg^^l^ have set out 
the formalism so we do not enter into it again here. 

The ground state that we access through MC is equiv- 
alent to the UHF result, but with no assumptions about 
translational symmetry. It is easy to see some of the 
qualitative effects of dynamical fluctuations in xrii and 
that we have neglected, at T = 0. These would 
(i) convert the U < Ud PM to a correlated metal, (ii) in- 
troduce quantum spin fluctuations in the large U AFI, 
and (iii) possibly shift Ud to a larger value (since the 
correlated metal competes better with the local moment 
phase). The intermediate window 'spin glass' that ap- 
pears within the static approximation might be converted 
to a spin liquid with slowly fluctuating moments. A re- 
cent calculation on the triangular lattice demonstrates 
how longer range and multi-spin interactions arise on a 
frustrated Mott insulator and can lead to a spin liquid 
ground statc'^'^. 

Our approach captures the correct thermal fluctua- 
tions of the nii, without any assumption about LRO in 
the background. This in turn allows us to capture a Tc 
that has the qualitatively correct U dependence. With 
growing temperature, but staying at T ^ [/, these clas- 
sical thermal fluctuations should reasonably describe the 
magnetic background, and its effect on the electrons. 

While our solution of the FCC Mott problem involves 
approximation, and real materials usually require inter- 
actions and degrees of freedom beyond the Hubbard de- 
scription, our results suggest the organization of a broad 
class of experiments, (i) We find that a non coplanar spin 
configuration dominates the Mott phase, and pressure in- 
duced metallisation leads to a weak moment 'spin frozen' 
state with short range non coplanar correlations. This 
is consistentP^I^ with observations on GaTa4Se8, FCC 
A3C60, and double perovskite materials, (i i) Beyond the 
pressure driven IMT the materials exhibitP^HMIli] ygj-y 
high p(0), and dp/dT < 0. Our results show how this 
can arise from the presence of disordered local moments 
strongly coupled to the itinerant electrons, leading to 
large scattering. In the Ga cluster materials these local 
moments emerge from correlation effects, while in the 
pyrochlored^ they are already present as localised aig 
electrons, (iii) The very recent optical measurement!^ in 
GaTa4Se8 shows pronounced non Drude character in the 
metal near the IMT. This is consistent with our optics 
results for the PM to PI transition. In fact we have a 
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quantitative description of the Mott transition in these 
material^^ within our present framework. 

The correspondence above allows us to make two 
concrete predictions: (a) The frustrated Mott sys- 
tems should have a wide pseudogap regime beyond the 
insulator-metal transition, persisting to T = 0. These 
should be visible in tunneling and photo-emission spec- 
tra, (b). The thermally induced shift of single particle 
spectral weight would be extremely asymmetric. Weight 
at low positive frequencies is shifted to the scale oiu ^ U, 
while the negative frequency spectrum remains almost 
unaffected. 

We have not probed the anomalous Hall response due 
to flux like correlations and plan to present them sepa- 
rately. Some of the FCC Mott materials have a super- 
conducting instability at very low temperature, ~ lOK. 



We have not treated that aspect, and need a method that 
retains the dynamics of the nii to do so. 

Conclusions: We have provided the first comprehen- 
sive study of the Mott transition on the geometrically 
frustrated FCC lattice. The magnetic frustration leads to 
a seemingly 'two fluid' state of itinerant electrons and dis- 
ordered local moments between the paramagnetic metal 
and the antiferromagnetic insulator. The disordered 
phase involves a large residual resistivity, non Drude op- 
tical response, and a single particle pseudogap. The tem- 
perature and interaction dependence that we uncover al- 
lows a common conceptual scheme for a wide variety of 
materials. 

We acknowledge use of the HPC clusters at HRI. PM 
acknowledges support from a DAE-SRC Outstanding Re- 
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Supplementary information: 

A. Magnetic phases on FCC lattice 

In the cubic lattice notation, the primitive lattice 
translation vectors for FCC are Ai = a(0, 1,1),A2 = 
a(l,0,l),A3 = a(l,l,0), so that the nearest neighbour 
distance is dfcc = V2a. All the points on FCC lattice are 
expressed in integer units of these, i.e., X = "i^i — 

a(n2 + 71-3, 71-3 + rii, rii + 722). Each site has 12 neighbours 
X + S, where 5={±1, ±1, 0), (±1, 0, ±1), (0, ±1, ±1). The 
'flux' phase mentioned in the text, is described by the 
following formula 



m(X) =TO(e'^i-^,e*'^ 



Where, m is magnitude of the vector, and Qi = (^, 0, 0), 
Q2 = (0, f , 0), Qs = (0, 0, ^). Its a non-coplanner phase. 
On the other hand, the 'C-type' phase is a coUinear one, 
described by the formula 

m{X) ^ m(0,0,e*'5<=-^) 

Where, Qc = (^,^,0). It consists of alternating ferro- 
magnetic lines. Both the 'flux' and 'C-type' phases are 
shown in Figj5]in the top and bottom panel respectively. 



B. Derivation of the effective Hamiltonian 

Consider the single band Hubbard model as starting 
point 

= Ho - fi^Ui + ni-fUii 

i i 

and implement a rotation invariant decoupling of the 
Hubbard term as follows. First, we write 




FIG. 5: Color online: The 'flux' (top) and 'C-type' (bottom) 
phases, drawn on cubic lattice motif. The green spheres are 
the FCC lattice points, and arrows indicate the direction of 
mi moments. The tetrahedra highlights the geometrical frus- 
tration. 



where — rii^ -I- n.4 is the charge density, Si — 
\ So p '^ia^apCip = 2(Ti is the local electron spin oper- 
ator, and rhi is an arbitrary unit vector. 

The partition function of the Hubbard model is 



Z = j D[c,c]e- 

S ^ drCir) 
Jo 



C = ^Qo-(r)9rQa(r) -f ff(T) 



We introduce two space-time varying auxiliary fields for a 
Hubbard-Stratonovich transformation: (i) ^^(t) coupling 
to charge density, and (ii) Ai(T)mi(r) = mi(T) coupling 



to electron spin density (A^ is real positive). This allows 
us to define an SU(2) invariant HS transformatiorpJ^, 



J 47r2[/ ^ 

The partition function now becomes: 

_ /■ T-T dcjdcidcjiidxai {- c{r)) 
^ y 11 47^2;/ ^ 

i 

C{t) ^ ^ Cia{T)drCia{T) + Hq{t) + /I^^f (0j (t) , Hlj (t) ) 
ia 

^int — ^ ^ 



— -f i(i>ini -I- - 2va.i ■ Si 



As mentioned in the text, to make progress we use two 
approximations: (i) neglect the time (r) dependence of 
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the HS fields, (ii) replace the field (f>i by its saddle point 
value {U /2){ni) = C//2, since the important low energy 
fluctuations arise from the 1x1^. Simplifying the action 
with above substitution, we get the effective Hamiltonian 

i i i 

where jl = /i— J7/2. For the sake of convenience we rescale 
— >■ ^nii, so that the nij is dimensionless. This leads 
to the effective Hamiltonian that we used in the text: 

Heff = Ho- fiYln,-^Y^ni,-a, + jJ2m^^ 

i i i 

The partition function can be written in terms of H^ff 

Z = J Vm.Tre^e^e-f^"^" 

For a given configuration {mi} the problem is quadratic 
in the fcrmions, while the configurations themselves are 
obtained by a MC sampling as discussed in the text. 



C. Optical conductivity 

The conductivity of is calculated as follows (ref!^, us- 
ing the Kubo formula: 

Jx = [^(c|.cr'^-'+a+t/,CT — he) 

i.cr 

Where, is the current operator and, the coefficient 
2 

A = ^ = y=^^ — . o'o=^ is the scale for conductivity 
with dimension of conductance. 

Ua — /(cq) is the Fermi function, and Cq. and \a) are re- 
spectively the single particle eigenvalues and eigenstates 
of Heff in a given background {nii}. N is the number of 
sites (N — 12^ in our results). The results we show in the 
text are averaged over equilibrium MC configurations. 

The d.c conductivity is the — > limit of the result 
above. The experimental value of dfcc for the GaTa4Se8 
cluster compound is ^ 4.3A. 
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